22, 2010 15:10 WSPC/INSTRUCTION FILE tavecchio 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



GAMMA-RAY EMISSION FROM AGNS (SPECIAL FOCUS ON BL 

LAC OBJECTS) 



F. TAVECCHIO 

INAF-OAB, Via Bianchi 46 
1-23807, Merate, LC, Italy 
fabrizio.tavecchio@brera. inaf. it 

Received Day Month Year 
Revised Day Month Year 
Communicated by Managing Editor 

Blazars, radio-loud active galactic nuclei with the relativistic jet closely aligned with the 
line of sight, dominate the extragalactic sky observed at gamma-ray energies, above 100 
MeV. We discuss some of the emission properties of these sources, focusing in particular 
on the "blazar sequence" and the interpretative models of the high energy emission of 
BL Lac objects. 
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1. Introduction 

The extragalactic gamma-ray sky at energies above 100 MeV is mainly populated 
by blazars, radio loud active galactic nuclei (AGNs) whose emission is dominated 
by the relativistically boosted non-thermal continuum produced in a relativistic jet 
(bulk Lorentz factor T = 10 — 20) pointing (typical angles less than 5 degrees) toward 
the Earth. Blazars are among the most powerful and violently variable sources in 
the Universe: their emission covers all the electromagnetic spectrum, from the radio 
to the gamma-ray band, and their bolometric output (often concentrated in the 
gamma-ray band) can exceeds L > 10 48 erg s _1 . The variability timescale can be as 
short as few minutedSHS] m t eres t for these sources is driven by the possibility 
to probe the innermost regions and the surroundings of the relativistic jets, close 
to the acceleration and collimation region (d ~ 10 2 — 10 3 gravitational radii). 

In the past decade EGRET onboard CGRO detected almost 70 blazars. Of the 
106 sources in the recently released listP (LBAS) of the AGNs detected with high 
significance by the Large Area Telescope (LAT) onboard the Fermi Gamma Ray 
Telescope^ in the first three months of operation (August-October 2008), 104 are 
blazars, the remaining two sources being the radiogalaxies Perseus A and Centaurus 
A. Note that one of the sources classified as a blazar is a radio loud Narrow Line 
Seyfert 1 galaxy, the first object of this class revealed at gamma rayi^rSl. The 
gamma-ray emission can extend up to the VHE band (E > 100 GeV), at which 
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the observations are performed by the Cherenkov telescopes on the ground. Also 
in this band blazars dominate the extragalactic sky. The VHE emission can be 
used to derive limits on the density of the IR-optical backgrounds exploiting the 
modification of the spectrum due to the absorption through the pair production 
procesJ^ - ^. 

Due to the limited space, in the following I will focus the attention on some 
specific issues: i) the blazar sequence and its physical meaning; ii) the high-energy 
view of the the BL Lac objects and some problems possibly requiring a revision of 
the simple one-zone emission models. 



2. SEDs and the blazar sequence 

The spectral energy distribution (SED) of blazars displays two broad peaks, inter- 
preted as due to synchrotron and inverse Compton (IC) mechanisms in the popular 
leptonic models^ - "^ (see, e.g., Ref . HTl fo r alternatives hadronic models). As shown 
almost ten years ago by Fossati et alM^ 1 , averaged SEDs of blazars in different bins 
of radio luminosity show a trend with the radio (or the bolometric) luminosity: 




Fig. 1. Energy index of the gamma-ray spectrum, en~, versus the gamma-ray luminosity for all 
the blazars detected with high significance by the LAT instrument onboard Fermi in its three 
months of observations. Open circles represent BL Lac objects, open squares flat spectrum radio 
quasars, filled symbols sources detected in the TeV band. Solid lines connect points corresponding 
to different emission states of the same source. A trend between the two quantities is clearly 
visible, together with a sharp division between the two classes of objects occurring at a luminosity 
L 7 ~ 10 47 erg s — 1 . These properties confirms the expectation of the so called "blazar sequence" 
(see text for details). 
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specifically, flat spectrum radio quasars (FSRQs), the most powerful blazars, show- 
ing intense broad emission lines and sometimes a thermal bump in their optical 
spectra, have both peaks located at low frequencies, submm-IR and MeV for the 
synchrotron and IC peak, respectively. On the contrary, in the less powerful BL 
Lac objects (characterized by very weak emission lines), the two bumps peak at 
UV- X-ray frequencies and above 100 GeV. The existence and the validity of such a 
trend (dubbed "blazar sequence") has been hotly debated in the pasfl^l. A possi- 
ble caveat concerns the use of the EGRET data to derive the averaged gamma-ray 
spectra, since the limited sensitivity likely introduces a bias toward flaring states. 
One argument of the critics is the claimed existence of "blue quasars" , FSRQs with 
peaks at high frequencies, violating the sequence. How ever, so far none of the blue 
quasars candidates has been proven to be a robust case^^l^. 

The three-month LAT data confirm with a larger statistic!^ one of the basic 
expectations of the blazar sequence, namely the fact that FSRQs should have (on 
average) soft GeV spectra (with energy index a 7 > 1) whereas BL Lac should show 
hard spectra (a 7 < 1). As shown in Fig|l]this prediction is confirmed by the new 
data, showing not only the expected dicotomy in the spectral slope, but also a well 
definite trend between the gamma-ray luminosity and the slope of the gamma-ray 
spectrum. The small "thickness" (or even the very existence) of the correlation is 
even surprising, considering the huge variability possible for both quantities in the 
same source. This is shown for some selected objects by the solid lines line connecting 
different emission states of the same source. In Ref. 1161 we proposed that the divide 
in the gamma-ray luminosity (a good proxy for the bolometric output of the jet) can 
be associated to a transition in the properties of the accretion flow on the central 
BH, transiting from a standard S hak ura-Su nyae v disk at high luminosities, to an 
inefficient accretion regime 

(ADApED, ADIOS22). The transition occurs at the level 
of an accretion rate around m ~ 10~ 2 in Eddington units, in agreement with the 
expectation from the theory!^ and the analogous scenario proposed to account for 
the FRI/FRII division among radiogalaxiea^l. 

The existence and the characteristics of the blazar sequence call for a robust 
and simple explanation. Interestingly, the spectral sequence is associated to a (even 
more robust) sequence in some of the fundamental quantities of the emission region, 
as modeled with one-zone leptonic models^} In particular, the energy of the elec- 
trons emitting at the peaks of the SEDs shows a strong correlation with the energy 
density (radiative plus magnetic) in the source. The most direct interpretation is 
that in sources with a large energy density the severe radiative cooling prevents the 
electrons to reach high energy, determining a low energy frequency of the peaks in 
the SEDs. In sources in which electrons are characterized small radiative losses the 
acceleration mechanism can instead bring electrons to higher energies, determining 
large peak frequencies. Electrons in the jet of FSRQs are commonly assumed to 
have a rather high cooling rate, since the environment external to the jet is ex- 
pected to be rather rich of soft photons, determining severe radiative losses through 
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Compton emission. BL Lac objects, instead, are characterized by an environment 
poor of radiation (and also possibly by a radiatively inefficient accretion flow) and 
thus electrons can be energized to very high energy, the cooling being dominated by 
synchrotron and synchrotron-self Compton losses. We refer to Ref. [22] for a deeper 
discussion. 



3. BL Lac objects at high energy: GeV and TeV emission 

As discussed above, jets in BL Lac objects are expected to propagate in a "clean" en- 
vironment and the emission we observed from them is likely due to synchrotron and 
SSC mechanisms. Being less prone to depend on the details specifying the location 
of the emission zone (e.g. L. Stawarz and L. Costamante, these proceedings) and 
the conditions of the surrounding environment, the emission from BL Lac objects 
seems to be the ideal tool to investigate the physics of the particle acceleration and 
cooling processes acting in jets. Moreover, BL Lac objects are interesting sources 
since the high energy component the SED is located at the crossroad between the 
band covered by LAT (E < 100 GeV) and by ground based Cherenkov telescopes 
(E > 50 GeV). It is thus possible to probe the entire high-energy emission, fully 
covering the bump before and after the IC peak (e.g. Refl2"3"f. 

To investigate the global properties of the BL Lac objects revealed at high and 
very high energies we considered all the BL Lac objects in the LB AS list, including 



1 





±111 1 1 111 1 111 1 1 11 1 1 |i 1 1 111 1 ii| 1 mi mi pin 111 
= o 



0„ 

□ oR a 

□ □ 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 u 
=1= o 



2-i 1- 

PQ 

» -2 I- 



o 



F°^Q EE o 8 § 

!□ a a a -EE- Sir * 



A 
A 



-3 
-4 



= O TeV/LAT A 
L * TeV no LAT J 

= □ LAT no TeV a: 

t I I I I : 



_ O 

□ o □ 



o 

A 



I I I I t 



2 3 4 5 6 7 20 30 40 50 
Log 7 b 6 



Fig. 2. Magnetic field versus the Lorentz factor of the electrons emitting at the SED peak (left) 
and the Doppler factor of the emitting region (right) derived for the BL Lacs using the one-zone 
SSC model. While most of the sources are characterized by magnetic fields around 0.1-1 G, Lorentz 
factors in the range 7^ = 10 3 — 10 4 and 15 < <5 < 30, there are some notable exceptions, requiring 
larger 8 and 7;, and smaller magnetic fields. Different symbols indicate the instrument (LAT or 
TeV telescopes) the detected the source. 
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also 12 sources detected at TeV energies. To this sample we add the remaining 12 
BL Lac detected at TeV energies (before the end of 2009 September) not detected 
by LAT. 

For all these sources we assembled the SED using the high energy data, historical 
data and, when available the (often nearly simultaneous) optical- UV and X-ray data 
from the Swift satellite. We reproduced the SED with a simple one-zone synchrotron- 
SSC model, deriving some of the basic physical quantities. For simplicity, the model 
assumes a spherical homogeneous region and the relativistic electrons are assumed to 
follow a broken power law distribution in energy with the break at the Lorentz factor 
7b. In Figj2]we report the values of the magnetic field B versus 7b (corresponding to 
the energy of the electrons emitting most of the power, close to the SED peaks) and 
the Doppler relativistic factor, i5. Most of the sources occupy a region specified by 
a magnetic field in the range B = 0.1 — 1 G, 7b = 10 3 — 10 5 , 6 = 15 — 30. However, 
there are some notable exceptions, requiring very small magnetic fields (B ~ 10~ 2 ) 
G, large 7b (10 6 - 10 7 ) and/or 8 > 40. 

At a closer look these "outliers" can be divided into two classes: i) sources 
showing a moderately hard LAT spectrum with the synchrotron peak in the optical- 
UV band; ii) TeV BL Lacs not detected by LAT characterized by extremely hard 
(de-absorbed) TeV spectra. 

3.1. Hard GeV spectra and low frequency synchrotron peaks 

The SED of one of the BL Lacs belonging to the group i), PKS 1717+177, is given 
in Fig|3] (adapted from Refl24|). Besides the LAT "butterfly", calculated with the 
spectral data of the LBAS catalogu we report the datapoints derived by analyzing 
the LAT data with the standard procedureQ The points confirm that the spectrum 
continues to be hard up to the highest energies. 

For these kind of sources the large separation between the synchrotron (optical- 
UV) and the SSC (E > 50 GeV) peaks in the SED inevitably requires large Doppler 
factors'^. In the case of PKS 1717+177 we show the model resulting by using a 
moderate value of the Doppler factor, 8 = 20 (solid line): the fit of the LAT data 
is clearly poor. Using a somewhat extreme value of 8 = 50 (dashed line) clearly we 
obtain a better agreement with the LAT data, although the LAT spectrum would 
be better reproduced with a peak beyond 50 GeV. With a word of caution due 
to the non simultaneity of the Swift and LAT data, we conclude that the simplest 
version of the SSC model would clearly be in trouble with these sources as far as 
we assume that 8 cannot exceed 50. 

A viable solution to alleviate the problem represented by the large Doppler 
factors is offered by the model assuming that the emission zone is not a single 
homogeneous region, but jet is structured, with a fast (F = 15 — 20) core (or spine) 
surrounded by a slower (r = 3 — 5) layer. These models were originally introduced to 



a http : //f ermi . gsfc .nasa. gov/ssc/dat a/analysis/ 
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Fig. 3. SED of PKS 1717+177. Optical-X-ray data are from Swift (January 2009). The "butterfly" 
reports the power-law fit of the first 3 month of LAT data. The corresponding data-points, have 
been obtained using the the standard analysis procedure for LAT data. The solid line is the result 
of the one-zone synchrotron-SSC model with a Doppler factor 6 = 20, the dashed line, instead, 
has been calculated assuming a somewhat extreme value of <5 = 50. 



account for the discrepancy between the large Doppler factors inferred in some of the 
TeV BL Lacs and the slow speeds measured in their jets at VLBI scald^ and have 
been subsequently adopted to model the high energy emission of radiogalaxieiP^. 
The key ingredient is the radiative coupling between the spine and the layer: both 
components see the emission from the other relativistically boosted (because of the 
relative speed). The IC emission is therefore amplified with respect to the SSC case. 
In particular, if the emission from the layer contains more soft photons than that 
of the spine, the corresponding IC-scattered photons can reach higher energy than 
the SSC photons, whose energy is strongly limited by the rapid decline of the Klein- 
Nishina cross-section. The emission from the spine could then show a high-energy 
tail from this kind of external Compton component, accounting for the observed 
LAT spectrum. A full discussion, including more detailed models is in preparation. 



3.2. Hard TeV spectra 

Among TeV BL Lacs there are a handful of peculiar sources showing very hard TeV 
spectra (when the observed spectrum is properly deabsorbed to take into account 
the absoprtion by photons of the cosmic IR-optical background) . The SED of one 
of the most extreme cases, 1ES 0229+200, is reported in FiggR 

In the framework of the one-zone SSC model such hard spectra are difficult to 
obtain, since the decreasing of the scattering efficiency in the Klein-Nishina regime 
leads to soft spectra. However, extremely hard spectra can be obtained^ if the 
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Fig. 4. SED of the BL Lac 1ES 0229+200, assembled with historical, Swift and HESS data. 
The optical-UV data track the emission from the host galaxy (dotted line). The solid line shows 
the result of the calculation with a synchrotron-SSC model assuming a simple truncated power 
law form for the electron energy distribution, with a large value for the minimum Lorentz factor, 
7min = 8-5 x 10 5 . The dashed lines report the LAT upper limits assuming a photon index of T = 2 
or 1.5. 

emitting relativistic electrons have a large value of the minimum electron energy 
7 m i n rac 2 : in this case the synchrotron emission, below the typical frequency of the 
low-energy electrons, z^min — 3 x 10 6 i?7^ lin (5 Hz, will display the typical F{v) oc v 1 / 3, 
spectrum. The same hard slope will be reflected by the SSC emission below the peak. 

In Fig|4]we report the optical-UV and X-ray data obtained by Swift, showing 
that a hard spectrum below the soft X-ray band is required by the data, supporting 
the possibility discussed above. Clearly, another strong prediction of this scenario is 
that the level of the emission in the GeV band should be rather low. Consistently, 
all the sources featuring such hard TeV spectra have not been detected in the first 
months of Fermi. In a recent paper^S, the LAT collaboration reports a deep (5.5 
months of observations) upper limit for 1ES 0229+200 (that we plot in FigJH dashed 
lines), fully consistent with the theoretical expectations. If our scenario is correct, 
these sources would be the most efficient cosmic accelerator known, able to push 
virtually all the electrons at very high energy. 

4. Concluding remarks 

1) The "blazar sequence" is confirmed and even strengthened by the data from LAT. 
I would like to remark, however, that deeper data will start to probe a regime of 
luminosities populated be a mix of sources. In particular, blazars with small black 
holes or with jets slightly misaligned will introduce some confusion, making the 
discussion more difficult. 

2) The one zone synchrotron-SSC model for BL Lacs gives satisfactorily result for 
most of the sources when we use almost simultaneous data. However, for a real test 
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of the model it is important to obtain really simultaneous data. 
3) The BL Lac mainly populating the LAT sky are, obviously, those with a promi- 
nent GeV emission. However, there is an interesting population of BL Lacs, char- 
acterized by very hard GeV-TeV spectra that are not easily detected by LAT. Any 
discussion on the high-energy emission of BL Lacs based only on LAT would exclude 
these sources that could be the most efficient machines accelerating the relativistic 
electrons. 
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